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In this article we use a lattice-Boltzmann simulation to examine the effects of shear 
flow on a equilibrium droplet in a phase separated binary mixture. We find that large 
drops break up as the shear is increased but small drops dissolve. We also show how the 
tip-streaming, observed for deformed drops, leads to a state of dynamic equilibrium. 

1. Introduction 

The study of phase transitions under shear has been of great interest in the recent 
pastB&tffl. To understand the basic phenomena underlying this complex process 
we focus our attention on the behavior of a single equilibrium droplet in a two- 
dimensional binary fluid under shear flow. Despite the simplicity of the model 
system it shows rich behavior, both droplet break-up and droplet dissolution. We 
also suggest a new explanation of tip streaming observed in our simulations. 



lis problem has 
Experimental, 
and theoret- 



ic 5 



Consider first an immiscible drop subjected to a shear flow, 
been studied extensively since the original experiments by Taylor 
theoretical and numerical results are available in three dimensions 
ical and numerical EtE£| results in two dimensions. These approaches consider 
drops with a singular interface and a conserved volume. The drops are deformed by 
the shear flow while maintaining their volume. If the shear rate exceeds a certain 
critical value, which depends on the volume of the drop, the drop will break up. 
Conversely, for a given shear rate, there exists a volume above which the drop is 
unstable. We shall denote this volume Vt>. 

For a miscible binary mixture a similar break-up of droplets is observed if the 
droplets are large. However there is now a second volume scale Vd which sets a 
lower limit to the drop size. Vd corresponds to the minimum size of a nucleation 
seed. The reason for the existence of Vd lies in the free energy balance between 
the favorable creation of separate phases in the supersaturated mixture and the 
unfavorable creation of the interface separating them. Note that, because when a 
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shear is applied a drop deforms and increases its surface length, Vd will depend on 
the shear rate. For shear rates with V> < Vd there are no stable drops in thesystem. 

To study these phenomena we use a model developed by Orlandini et. al&3 for the 
isothermal flow of a binary mixture. In this approach a free energy which describes 
a binary fluid is chosen. The pressure tensor and chemical potential calculated from 
this free energy are then included in a lattice Boltzmann scheme for modeling fluid 
flow. The fluid obeys the Navier-Stokes and convection-diffusion equations and 
comes to an equilibrium corresponding to the minimum of the input free energy. 
The coexistence curve is correct and interfaces are extended in space as predicted 
by the Cahn-Hilliard theory. One advantage of this approach is that the free energy 
and the chemical potential can easily be calculated as local functions and compared 
to theoretical predictions. 

In the next section of the paper we outline the thermodynamics of the binary 
fluid and describe the extensions to the lattice Boltzmann approach needed to treat 
shear flow. Results for the break-up of a large droplet are presented in section 
three. In section four we obtain an estimate for the volume Vd below which small 
droplets dissolve and discuss the effect of shear flow on the dissolution. Section five 
discusses tip streaming, the loss of material from the tips of the deformed droplet, 
and summarizes the results of the paper. 

2. Method 

We simulate a binary fluid comprising two components A and B, say. A-A and 
B-B interactions are zero but there is an A-B repulsion Xuatib where ua and ub 
are the number densities of A- and B-particles respectively. This system can be 
described by the Landau free energy functional 



* = J dr(^,n,T) + -(V^) 2 ) (1) 

where T is the temperature, n = ha + ns, <p = tia — ns and n is a measure for the 
interface free energy (surface tension). The free energy density of the homogeneous 
system is 

iP(<p,n,T) = ^(l-^-Tn 

+ |(n + ,)ln(^)4(n-,)ln(^). (2) 

The Navier-Stokes and convection-diffusion equations for the fl uid are simulated 
using a lattice Boltzmann approach described in detail in IIjEJ. Here we restrict 
ourselves to a description of the way in which shear flow was implemented. 
We consider a linear shear flow with velocity 



u x \ _ ( Gy 
u„ J ~ I 



(3) 
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Figure 1: Sketch of a deformed drop in a simple shear flow. For small shear rates 
the drop has the form of an ellipse with axes 2 A and 2B. The ellipse is inclined to 
the direction of the shear flow by a shear strength dependent angle a. 

where G is the shear rate. If the fluid is homogeneous it is expected that eqn.(|3|) 
describes the velocity field of the whole fluid. Inserting a droplet will disturb the 
velocity field locally but eqn.(^) gives the far field solution. 

To simulate the shear it is necessary to introduce boundary conditions that force 
the flow. After each streaming step we replace the collision step at the boundary 
by a step that defines the variables of the lattice Boltzmann scheme to take values 
that correspond to the required value of the velocity and densities. 

Two different kinds of boundary conditions have been implemented: 

periodic: The top and bottom edges of the lattice at y — iyb are the boundaries 
and the velocity is constrained to be u = (Gy, 0). The side boundaries have 
periodic boundary conditions. This corresponds to a shear confined by two 
moving walls acting on a periodic array of drops. 

forced: All the edges of the lattice are forced to have u = (Gy, 0). This eliminates 
the effects of periodic images of the drop. The local values of n and (f are 
replaced by their mean value averaged over the boundary. These boundary 
conditions generalize for more complicated forced flows, for example hyper- 
bolic shear flow. 

For a homogeneous system both boundary conditions lead to the velocity profile 
of eqn.(^|) to within machine accuracy. In the presence of a drop both boundary 
conditions are expected to give the same results for an infinite lattice. A comparison 
of the two different boundary conditions therefore gives a measure of the effect of 
the periodic images on the drop. 

If a drop is placed in a shear flow it will be deformed by the forces acting on it. 
The drop elongates and turns to lie at an angle a to the flow until in the steady 
state the restoring force due to the surface tension balances the shear forces acting 
upon the drop. This situation is sketched in figure [l| For small deformations the 
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Figure 2: a) Deformation of a drop D and b) the tilting angle a for periodic ( — ) 

and forced ( ) boundary conditions against the shear rate G. The undeformed 

drop has a radius of 8 lattice spacings and the lattice size is 60x30. 



drop approximates well to an ellipse and its deformation can be defined as 

D 



A + B 



(4) 



where A and B are the major and minor axes respectively. For drops of constant 
volume the deformation and inclination angle depend only on a dimensionless quan- 
tity, the capillary number§ 

Ca = ^ (5) 
a 

where v is the viscosity, a the undeformed drop radius and a the surface tension. 
Throughout this paper we take v = 1/6 and a — 0.046 (n — 0.002). This corre- 
sponds to an interface width « 3 lattice spacings. 



3. Break-up 

Typical results for large drops (V » Vd) are shown in figure || where D and a 
are plotted as function of the shear rate. Results for forced and periodic boundary 
conditions are compared. The results presented are for a lattice of size 60x30 with a 
drop of initial radius 8. They were obtained by equilibrating the fluid at each data 
point and then increasing the shear, re-equilibrating to give the next point, and so 
on. 

There is a linear dependence of the deformation on the shear rate for small shear 
rates, followed by a more rapid deformation as the shear increases, and finally break 
up. The different boundary conditions lead only to small quantitative differences in 
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Figure 3: a) An equilibrium drop at G — 0.002246, b) the same drop at G — 
0.002251 in the first stage of breaking up. The drop in b) breaks up into two drops 
of equal size at larger times. The arrows show the velocity at every third lattice 
point in each direction. 

the result. The results are in qualitative agreement with the results expected from 
comparison with three-dimensional experiments. For a drop of constant volume 
and given ratio of the viscosity of the drop^nd the surrounding fluid the deforma- 
tion depends only on the capillary numberu. The drops studied here do not have 
a constant volume and therefore deformation at break up should depend on the 
parameters of the system but it is reasonable to expect only a weak dependence. 
Indeed we find D\, ps 0.65 for drops with undeformed radii 8,13,20 in those cases 
where they break up. 

It is known! that the second curvature is very important for the rupture of 
three-dimensional drops. This mechanism does not exist in two-dimensional sys- 
tems. Therefore we felt it was important to check the existence of the break up 
carefully, particularly as it is well knownB that a sudden change in shear strength 
can lead to rupture long before the critical shear rate. We therefore performed care- 
ful numerical simulations where we saved a stationary solution and then increased 
the shear flow. If the drop ruptured instead of reaching a stationary state we loaded 
the old configuration and increased the shear rate by only half the previous amount. 
We iterated this step until the increase in shear rate was smaller than a lower bound. 
The shear rate never grew larger than a previously rejected shear rate showing that 
the break-up is not due to non-equilibrium effects. We then decreased the shear 
rate to check that the system is in equilibrium at every point. No hysteresis effects 
were observed showing that this is indeed the case. 

It should be pointed out, however, that the rupture mechanism is very different 
in two and three dimensions. Three dimensional drops break up into two main 
drops and satellite dropsB. We observe that two dimensional drops break up into 
two drops of equal size as shown in figure |^. 

4. Dissolution 

For given ua, rig, T, if a binary fluid lies in the coexistence region it will separate 
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Figure 4: The free energy plotted as function of the drop radius for a) a lattice 
with 10000 points for concentrations that without surface effects would correspond 
to drops of radii R° =15,16,17,18 and 19; b) a system with R° — 19 and lattice 
sizes 10000, 12000, 14000, 16000, 18000 and 20000; c) a system of size 10000 and 
R Q = 19 for shear rates G =0, 0.001, 0.002, 0.003, 0.004 and 0.005. 

into two phases with concentration differences ipi and ip2, say. Each phase occupies 
a fraction of the total volume of the system which is determined by the densities 
and temperature. For a finite system the interface between the two phases results 
in a finite positive contribution to the free energy and, if this is too large relative 
to the gain in free energy due to the phase separation, the drop is unstable. 

For the model simulated here it is possible to obtain an estimate of the droplet 
volume below which dissolution will occur. Ignoring the interface curvature the 
surface tension for an interface orthogonal to the z-direction iaij 




Hence, using (|l]) the free energy of a drop of radius R in a volume V is 

F = nR 2 ^(ip 1: n, T) + (V — wR 2 )i}(<p2, n, T) + 2irRcr. (7) 

This should be minimized with respect to R, ipi and ip2 where one variable can be 
eliminated by the constraint 

irR 2 cpi + (V - TrR 2 )ip 2 = (p. (8) 

In figure ||a the free energy, relative to that of a homogeneous system, is plotted 
as a function of the drop radius for different concentrations of A and B particles. 
The concentrations are chosen such that without the surface effects drops of radii 
15,16,17,18, and 19 would minimize the free energy on a lattice of size 100x100. We 
will denote these radii by R°. Without the surface free energy term there is only one 
minimum and this will always be reached. With surface effects the homogeneous 
phase (R = 0) is always stable and a finite deviation from it is needed to reach the 
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Figure 5: Deformation as function of shear rate for a drop of radius 8 for systems of 
size a) 60x30, b) 70x40. The drop in the smaller system breaks up while the drop 
in the larger system dissolves. Because the mass of the drop decreases strongly 
at high shear rates its deformation decreases under the increase of shear before it 
completely dissolves. 

global minimum. This corresponds to the metastability of some regions of phase 
space where a finite nucleation barrier prevents immediate phase separation. 

The lowest free energy curve in figure ||a corresponds to R° = 19. The figure 
shows that the effect of the surface free energy is to shift the minimum from 19 to 
about 16.7. For smaller R the minimum becomes a local minimum at R ~ 18. In 
systems with fluctuations this will eventually lead to the drops dissolving. However 
the lattice Boltzmann simulations reported here do not include noise and a drop in 
a local minimum will be stable. For concentration ratios that lead to a graph that 
has no local minimum (R° ~ 15.5) the drop will dissolve. However we observe that 
the initial dynamics are very slow. 

The radius below which drops dissolve depends not only on the size of the 
drop but also on the total volume of the system. Results for the concentration 
corresponding to R° — 19 for different system sizes are shown in figure ^b. Droplets 
are less stable in a larger system because more material from the drop is needed to 
change the concentration outside the drop. 

When shear is applied the droplets deform. They have a larger interface so that 
the surface contributions are increased. For small shear rates the dependence of the 
deformation on the shear rate is approximately linear and we can use 



where R is the radius of the undeformed drop with the same volume. Approximating 
the shape of the deformed drop by an ellipse the length of the interface is given by 



D 



10GR 



(9) 




(10) 
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Figure 6: a) The density difference ip and b) the chemical potential /i for a highly 
deformed drop. 

where E(k) denotes the complete elliptic integral of the second kind. This result 
can be used to estimate the effect of shear flow on the free energy. The results for 
R° = 19 for different shear rates are shown in figure ^c. The minimum in the free 
energy vanishes for high shear rates and no stable drops can exist. This corresponds 
to the dissolution of a drop under shear. 

The effects predicted by the thermodynamic theory are observed in the simula- 
tions. In figure |5| results for the simulation of a drop of R° = 8 on lattices of size 
60x30 and 70x20 are shown. For the smaller lattice Vb < Vd and the drop breaks 
up. For the larger lattice the drop can exist to larger values of the shear, but loses 
mass. Finally as the shear is increased further it dissolves. 

It is interesting to note that the theory predicts that every drop will dissolve 
under shear in a sufficiently large system. That is because the mass loss from 
the drop has to change the concentration outside the drop in order to reach a 
new equilibrium. This can, however, be very difficult to observe because there is 
a separation of time scales. The time tp> needed for diffusion to equilibrate the 
system scales as tn ~ L 2 where L is the length of the system. This result is altered 
under shearE3 and in the limit of very large shear rates tp, ~ L* . The time tp for 
the system to reach a steady flow and a new deformation of the bubble scales as 
tp ~ L. Therefore in large systems the time for the deformation of the drop will be 
fast compared to the time for it to dissolve. 
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5. Discussion 

In our previous discussion we have assumed that thermodynamic arguments for 
equilibrium can be carried over to the stationary states of dynamic systems. This 
is, however, not necessarily the case. We observe a deviation from this assumption 
in a phenomenon which we will call "tip-streaming" . Mass is pulled from the ends 
of the drop by the shear flow. The depletion of the concentration in the drop leads 
to a reduction of the chemical potential in the drop. The non-constant chemical 
potential leads to a diffusion current acting in the direction of the chemical potential 
gradient which returns material to the drop. Equilibrium is reached when the dif- 
fusion into the drop balances the tip streaming. The chemical potential in dynamic 
equilibrium is shown in fig.|^b. Note that this is a state of dynamic equilibrium not 
thermodynamic equilibrium as fi is not constant. 

For a small diffusion constant a large chemical potential difference is needed 
for equilibrium and hence a large amount of mass is pulled from the drop. This 
mechanism for the dissolving of a drop is distinct from the free energy driven mech- 
anisms explained earlier in the paper. It depends on the diffusion constant which 
is a dynamical quantity that does not enter the free energy. For infinite diffusion, 
however, material is immediately returned to the drop and no difference in the 
chemical potential is set up. A similar tip=-streaming was observed in experiments 
but was interpreted as a surfactant effect 13. 

Halliday et.al. performed simulations on a droplet in a binary fluid using a 
derivative of the the GunstensenEH algorithm. They obtain a similar break-up be- 
havior but there are some important differences in the results. In particular in the 
simulations r epo rted here all equilibrium drops are convex (see figure |^) whereas 
Halliday et.altn report stationary drops that are constricted in the middle. They 
also observe smaller inclination angles a at break up. The discrepancies warrant 
further investigation. Another feature not reported by Halliday et.al. is the disso- 
lution of drops under shear. This may be due the fact that they work in a region 
where the two components are more strongly separated. Nevertheless small droplets 
should still dissolve. We caution that the time-scale for dissolution can be much 
larger than that for the deformation of the droplet to reach an equilibrium value. 

Goldburg and Mini performed experiments on nuclcation in a binary mixture 
in the presence of shear. They observed the vanishing of drops under the influence 
of shear as a sharp transition. This transition may be interpreted as corresponding 
to the point where VJ, = Va- 

To conclude, in this paper we have shown that drops under increasing shear flow 
can either break up or dissolve. We explained this behavior with thermodynamic 
arguments. Tip-streaming was shown to lead to a state of dynamic equilibrium for 
deformed droplets. 
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